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We present near-infrared imaging of a sample of 57 relatively large, Northern spiral 
galaxies with low inclination. After describing the selection criteria and some of the 
basic properties of the sample, we give a detailed description of the data collection and 
reduction procedures. The K s X — 2.2/im images cover most of the disk for all galaxies, 
with a field of view of at least 4.2 arcmin. The spatial resolution is better than an 
arcsec for most images. We fit bulge and exponential disk components to radial profiles 
of the light distribution. We then derive the basic parameters of these components, as 
well as the bulge/disk ratio, and explore correlations of these parameters with several 
galaxy parameters. 
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1 INTRODUCTION 



Near-infrared (NIR) imaging of galaxies is a better tracer 
of the underlying stellar structure than optical imaging, for 
two reasons. Firstly, the effects of extinction by dust in the 
galaxy are much reduced (an order of magnitude from V to 
K), and secondly, the NIR emission is a better tracer of the 
older stellar populations, which contain most of the stellar 
mass. The combined effect of these two factors is that NIR 
emission is a more accurate tracer of the mass in galaxies. 
This is also why the fraction of galactic bars, mostly de- 
void of current star formation (SF), is higher as determined 
from NIR than from optical surveys (~70% versus ~ 55%, 
e.g., Sellwood & Wilkinson 1993; Mulchaey & Regan 1997; 
Knapen, Shlosman & Peletier 2000; Eskridge et al. 2000). 

The main technical problems in NIR imaging of spi- 
ral disks are the high background at 2.2 micron, and the 
relatively small size of NIR array detectors. Detectors with 
1024x1024 pixels, such as the ones used for the present pa- 
per, have only been in use since about 1996, and many of 
the cameras built around such arrays have been constructed 
to exploit the higher spatial resolution in the NIR, and/or 
the possibilities of adaptive optics (e.g., KIR on the CFHT, 
Doyon et al. 1998; UFTI on UKIRT, Roche et al. 2002). So 
whereas optical surveys of galaxy disks have been performed 
for decades, resulting, among many other results, in galaxy 
catalogues such as the RC3 (de Vaucouleurs et al. 1991), NIR 
surveys covering the whole disk of nearby spiral galaxies are 



only now starting to be published (e.g., 2MASS: Skrutskie et 
al. 1997, Jarrett et al. 2003; Seigar & James 1998a, 1998b; 
Moriondo et al. 1999; Mollenhoff & Heidt 2001; Eskridge 
et al. 2002; this paper). The NIR cameras we used for the 
presently described project have been designed as a com- 
promise between sampling the best seeing conditions and a 
relatively wide field of view (INGRID on the 4.2 m William 
Herschel Telescope, Packham et al. 2003), or specifically for 
wide- field imaging (PISCES on the 2.3 m Bok telescope, Mc- 
Carthy et al. 2001). The presently described survey offers 
several advantages over most of the other studies mentioned 
above, namely by offering a set of high-resolution (the me- 
dian value for the FWHM seeing across our sample is 0.90 
arcsec), well-sampled (pixel size 0.24 arcsec for all but 4 im- 
ages), reasonably deep images of a representative sample of 
nearby, not highly inclined, Northern spiral galaxies. 

Decomposition of galaxy light in bulge and disk com- 
ponents has received renewed attention in recent years (e.g., 
Andredakis, Peletier & Balcells 1995; Byun & Freeman 1995; 
de Jong 1996a; Graham 2001; Mollenhoff & Heidt 2001). 
One of the main reasons for this renewed attention stems 
from the realisation that not all spheroids have R* de Vau- 
couleurs (1948) type profiles, but often are better fitted by 
R~ Sersic (1968) type profiles. Furthermore, it has been sug- 
gested that the shape of the bulge and its relation to the disk 
may hold clues about the formation history of the galaxy 
components (e.g., Courteau, de Jong & Broeils 1996; An- 
dredakis 1998; Graham 2001; Trujillo et al. 2002). It is now 
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generally believed that early-type spiral galaxies have bulges 
with de Vaucouleurs-like profiles (i.e. Sersic n values of order 
4) and late-type spiral galaxies bulges with exponential-like 
profiles with n ~1, probably with a continuum in between. It 
has been suggested that the exponential bulges in late-type 
spiral galaxies are the result of bar instabilities driving gas 
toward the centre and scattering star in the vertical direc- 
tion by buckling resonances, while the de Vaucouleurs type 
bulges are the result of accretion and merging events, prob- 
ably early in the history of the galaxy (e.g., Courteau et al. 
1996; Aguerri, Balcells & Peletier 2001). An alternative view 
was put forward by Andredakis et al. (1995), who proposed 
that the outer parts of bulges were strongly affected by 
the disk in late-type spiral galaxies, resulting in bulges with 
lower n values. Andredakis (1998) found that by growing a 
disk adiabatically on an existing bulge in a collisionless 
N-body simulation, he could indeed lower the n value of the 
resulting bulge, but not to values lower than n = 2. 

In this paper, we present NIR imaging in the 2.2 micron 
K s band of a complete sample of 57 nearby, not highly in- 
clined, Northern spiral galaxies of all types. We also present 
fits to the separate bulge and disk components on the basis 
of radial surface brightness profiles. Optical broad (B, I) and 
narrow (Ha) band imaging of all galaxies in our sample will 
be presented in Paper II (Knapen et al. 2003). In Paper III 
(J. H. Knapen, in preparation) we study the morphology of 
the Ha emission in the central few kiloparsecs of the sample 
galaxies to explore the relative frequencies of, for instance, 
central peaks or nuclear rings among the sample. Correla- 
tions with morphological type or nuclear non-stellar activity 
will then be explored. In further papers, we will use ellipse 
fits to the NIR images to explore the parameters of the bars 
in the sample galaxies, and will use disk scale lengths to in- 
vestigate whether, and in what fashion, these are affected by 
emission from young stars in areas of current star formation. 
The currently presented data set offers many more oppor- 
tunities for study, and as examples we mention the study 
by Block et al. (2001) who derived gravitational torques due 
to bars from our NIR images, and the one by Jogee et al. 
(2002), who used one of our NIR images, of NGC 5248, to 
argue for the dynamical connection of the different compo- 
nents identified in that galaxy. 



2 THE SAMPLE 
2.1 Sample selection 

Our sample of galaxies was selected from the list of 
Elmegreen and Elmegreen (1987), which contains 708 galax- 
ies from the Second Reference Catalogue of Bright Galax- 
ies (RC2; de Vaucouleurs et al. 1976). The Elmegreen and 
Elmegreen sample was selected to included all those galax- 
ies listed in the RC2 with declination 8 > —35°, inclina- 
tion i < 60° and inclination-corrected diameter at 25 mag 
arcsec -2 , D25, greater than 2 arcmin. Our sample of 57 
galaxies was extracted from this large sample to give a 
smaller subset, primarily for analysis of the spiral arm prop- 
erties, as follows: 

• All galaxies of 4.2 arcmin and above in diameter (-D25) 
were selected to ensure that we could isolate the spiral arms, 
even at lower spatial resolutions. 



• Galaxies with 5 < —20° were excluded, ensuring all 
galaxies could be observed from the Northern hemisphere. 

• Galaxies of inclination higher than 50° , where we would 
have difficulty isolating the spiral structure, were excluded. 

This gave a sample of 57 galaxies (Table 1) of which six 
have a diameter larger than 10 arcmin. The limitations of 
galaxy catalogues are various, due to their adopted criteria 
of inclusion and their degree of completeness. The RC2 from 
which the original sample was chosen is, as implied by its 
title, magnitude limited. No attempt has been made by us to 
investigate selection effects in the sample, but we note that 
the range of galaxies included in our sample is by nature 
limited in apparent diameter, brightness, and thus also in 
redshift z. This is a factor which must be taken into account 
with any conclusions drawn from the data. 

Our sample of spiral galaxies covers the full range of 
morphological types from SA to SB, the full range of Hub- 
ble types, and the complete range of arm class from 1 to 
12 (Elmegreen and Elmegreen 1987)*. The sample contains 
normal and active galaxies; isolated, companion and clus- 
ter member galaxies; disturbed and undisturbed galaxies; 
barred and un-barred galaxies; galaxies with and without 
extensive massive star forming regions; and galaxies with 
and without nuclear rings or other nuclear structure. Some 
of the most important of such basic parameters of the sam- 
ple galaxies are given in Table 1. Stedman & Knapen (2001) 
present graphically the distribution of the sample galaxies 
as a function of morphological type, arm class, bar type, 
disk diameter, ellipticity, and systemic velocity, and we thus 
refer to the figures in that paper. 

Of the sample galaxies, 23% are barred (RC3 type SB), 
26% are un-barred (SA) and 51% are mixed (SX), giving 
a total of 74% of galaxies with bar structure. This is con- 
sistent with the findings of, e.g., Knapen et al. (2000), who 
estimate a bar fraction of approximately 2/3rds for all galax- 
ies. Again based on classifications from the RC3, 25% of the 
sample galaxies possess an inner ring structure, 30% pos- 
sess an inner s-shaped structure and 45% are mixed. The 
sample covers the full range of Hubble types from SO/a to 
Sm, with a marked excess of type be and c galaxies. In 
fact, these two classes of open armed spirals contribute 40% 
of the total. Every arm class is represented in the sample, 
with 45% of the galaxies in classes 7, 8, 9 and 12, denoting 
grand design spiral arm structure. Our sample galaxies are 
fairly evenly distributed in recession velocity, lying mostly 
between 750 kms _1 and 1750 kms -1 , with a few galaxies at 
velocities between 2000 and 3000 kms -1 . We checked the 
NASA/IPAC Extragalactic Database (NED) for nuclear ac- 
tivity in our galaxies. As listed in Table 1, 27 galaxies (52%) 



* Elmegreen & Elmegreen (1987) define arm classes as follows: 
arm class l=chaotic, fragmented, unsymmctrical arms, 2=frag- 
mcntcd spiral arm pieces with no regular pattern, 3=fragmentcd 
arms uniformly distributed around the galactic centre, 4=only 
one prominent arm; otherwise fragmented arms, 5=two symmet- 
ric, short arms in the inner regions; irregular outer arms, 6=two 
symmetric inner arms; feathery ringlike outer structure, 7=two 
symmetric, long outer arms; feathery or irregular inner arms, 
8=tightly wrapped ringlike arms, 9=two symmetric inner arms; 
multiple long and continuous outer arms, and 12=two long sym- 
metric arms dominating the optical disk 
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are listed as showing some type of nuclear activity, of which 
almost half (13) are Seyferts and most others LINERs. 

We conclude that although the sample is limited by its 
selection criteria, most importantly a lower limit on diameter 
and is therefore not "complete" in a strictly statistical sense, 
it is internally consistent and a fair, representative sample 
of all galaxies within the selection parameters, i.e., bright, 
relatively face-on spiral galaxies. 



3 OBSERVATIONS AND DATA REDUCTION 
3.1 INGRID observations 

Images of all but four of our sample galaxies were obtained 
through the K s filter with the INGRID NIR camera (Pack- 
ham et al. 2003) at the Cassegrain focus of the 4.2 m William 
Herschel Telescope (WHT). The INGRID camera is built 
around a 1024x1024 pixel Hawaii detector array (HgCaTe), 

and gives a projected pixel size of 0.242 arcsec^ with a field 
of view of 4.2 x 4.2 arcmin. Observations for this programme 
were made during a total of six allocated observing nights 
(May 14-16, 2000, Nov. 8 & 12, 2000 and Dec. 7, 2000), while 
additional images were collected during a number of service 
nights (Nov. 3, Dec. 4, 2000; March 9, Oct. 4, 2001; and 
Jan. 4, 2002). The image of NGC 3351 was obtained dur- 
ing the commissioning of INGRID, on the night of March 
18, 2000. Conditions during most of the nights were good, 
photometric and with good seeing. Whereas for some galax- 
ies images could only be obtained on one night, many were 
observed during different nights. The images of worst qual- 
ity in terms of sky subtraction are of a number of galaxies 
around RA=14h, which could only be observed during our 
May 2000 nights, when the observations were hampered by 
clouds. Total exposure times range from 6 to 100 minutes on 
source, as listed in Table 2. Typical la background noise is 
20.5 mag arcsec -2 , depending on the observing time, down 
to 21.7 magarcsec -2 for our longest exposure, of NGC 7741. 
This compares favourably with the value given for 2MASS 
images of bright galaxies, of 20.0 magarcsec -2 (Jarrett et 
al. 2003). 

The observing strategy consisted of small blocks of 
galaxy observations, each of some 5 minutes duration, in- 
terspersed with similar blocks of exposures of a nearby area 
of background sky, for an equal amount of time. The blocks 
consisted of 4 pointings, where in each pointing the galaxy 
centre was displaced by about 10 arcsec on the array. Each 
of these pointings in turn was built up from 3 to 5 individ- 
ual exposures of 20 to 12 seconds, co-averaged to produce 
one FITS file after 60 second total exposure time. The high 
background in the K s filter, variable and mostly dependent 
on the temperature of the outside air and the telescope, set 
the limit to the length of the individual exposures. The so- 
called dither pattern of 4 pointings serves to eliminate the 
bad pixels during the reduction, while the background sky 
exposures are paramount for sky subtraction (see below). 

Whenever conditions were photometric, standard stars 

t The pixel size of INGRID is 0.238 arcsec as of March 22nd, 
2001. The few images we obtained after this date were re-griddcd 
to the original pixel size of 0.242 arcsec. 



from the list of Hawarden et al. (2001) were observed at 
regular intervals during the night. 

3.2 INGRID data reduction 

The data reduction was performed using a combination of 
standard and purpose-built IRAF and IDL tasks. The main 
reduction script has been adapted from a set of IDL proce- 
dures originally written by F. Rigaut and R. Doyon to re- 
duce images taken with the KIR camera on the CFHT. The 
first step is the creation of a background frame from offset 
sky exposures taken before and after the galaxy frames. The 
different background frames are combined using an iterative 
median combination algorithm in order to remove stellar 
images and bad pixels which may be present in the individ- 
ual exposures. In each of the (normally three) iterations the 
individual frames are compared with the median-combined 
frame produced from them, and any pixels which lead to 
peaks or troughs in the combined image of more than three 
times the noise are set to undefined in the offending single 
frame in the next iteration. 

The background frame is then subtracted from each 
individual galaxy frame, and the background-subtracted 
galaxy frames are median-combined into the reduced image. 
In this latter step, bad pixels are masked out, and flatfield- 
ing is done on the individual frames. Since the galaxy is 
at a slightly different position in each individual frame, the 
position of either the nucleus of the galaxy or of a star in 
the field is determined and used to shift the images before 
combining them. 

Where several images of the same galaxy were available, 
either taken during the same night or on different nights, 
these images were averaged after registering them using the 
positions of a number of foreground stars, and/or the nu- 
cleus of the galaxy. The procedures used for the alignment 
were very similar to those described in §3.4, below, for the 
registering of NIR and optical images. What we will refer to 
as "final" images are those which have been registered with 
the optical images of our data set (as described in detail in 
Paper II). 

Photometric calibration was performed by comparison 
of the final, combined, images with the reduced images as 
obtained on individual photometric nights. For those nights, 
calibration constants were derived from the observations of 
photometric standard stars. We measured the total flux in 
the two images in a circular aperture centred on the nucleus 
of the galaxy, making sure that exactly the same area was 
taken for both the photometrically calibrated and the final 
image. Correcting for background and pixel size, which are 
in some cases different in the photometric and final images, 
we thus calibrate the final images. 

3.3 PISCES observations and data reduction 

Images of four of our sample galaxies were obtained using 
the PISCES camera on the Bok 2.3 m telescope of the Stew- 
ard Observatory, on Oct. 17 (NGC 6946) and 18 (NGC 628, 
1068, 6184), 1999. The PISCES camera (McCarthy et al. 
2001) uses the same kind of array as INGRID, but gives 
pixels of (f.'5 on the sky. The details of the observations, 
data reduction and calibration procedures for these images 
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Table 1. Global parameters of the galaxies in the observed sample. Except where otherwise mentioned, all figures are from the RC3. 
Tabulated are NGC (column 1) and Messier numbers (column 2) of the sample galaxies; morphological type from the RC3 (column 3) 
and from the Carnegie Atlas of Galaxies (Sandage & Bedke 1994; column 4); nuclear activity from the NED (column 5); apparent major 
isophotal diameter measured at or reduced to surface brightness level /ig = 25.0-B-mag arcsec (D25; column 6); Elmegreen arm class 
(from Elmegreen and Elmegreen 1987; column 7), where the arm class ranges from 1 (flocculcnt) to 12 (grand design); inclination (i; 
column 8), from the ratio of the major to the minor isophotal diameter; position angle PA of the major axis of the disk (column 9); total 
magnitude Mg (column 10); mean heliocentric radial velocity as derived from neutral hydrogen observations (v; column 11); distance 
from the Nearby Galaxies Catalog (Tully 1988; column 12) and image scale (column 13), as derived from the distance. 
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are described in McCarthy et al. (2001) and only the es- 
sentials are repeated here. The PISCES camera has an 8.5 
arcmin circular field of view and the galaxies were mosaiced 
with a few central images and many overlapping images in 
the outer parts, interspaced with sky offset images. Flat- 
fields were created by differencing a number of early and 
late twilight images. The fiatfielded object images were sky 
subtracted by weight- averaged and clipped offset sky frames. 
The object images were corrected for image distortions and 
then optimally combined using cross-correlations for posi- 
tional offsets and using airmass coefficient calibrations for 
intensity matching. Effective exposure times range typically 
from about 1000 s in the outer regions to 2000 s in the cen- 
tre. Seeing on the combined images is 2.5 arcsec for the Oct. 
17 observation and 1.5 arcsec for those of Oct. 18. Calibra- 
tion was performed by observing standard star fields from 
the list of Hunt et al. (1998) at least 4 times a night at sev- 
eral airmasses, to allow accurate calculation of airmass and 
colour terms. 

3.4 Final data sets 

Since we are interested in comparing the NIR images with 
the optical images we have obtained for our sample galax- 
ies, we have used a special procedure to make the NIR im- 
ages directly comparable to the others. This procedure is 
described in detail in Paper II, where we also describe the 
optical observations and their reduction, but in outline the 
following was done. We determined the positions of three 
foreground stars (or failing that, two stars plus the galaxy 
nucleus) in each of the images to be combined, as well as in 
an astrometrically calibrated image which we obtained from 
the Digitized Sky Survey (DSS). From these positions, an 
IRAF script we developed calculates the translation, rota- 
tion, and scaling needed for each image. It also determines 
which image has the highest resolution (smallest pixels) and 
the largest extent. Finally, the rotation needed to place the 
image set in an RA, dec orientation is determined from the 
position of the stars in the DSS image. Images are then ro- 
tated, translated and scaled in such a way that the resulting 
images are correctly oriented, have the pixel size of the image 
with the smallest pixels, and the extend of that covering the 
largest area. Those images covering a smaller area are filled 
with zero-valued pixels. Even though the resulting images 
may contain many more pixels, and are thus much larger, 
this is the only way to ensure that no information is lost. 
Table 2 gives the size of the fully reduced K B images, but be- 
fore the operations described immediately above. The area 
given in the Table thus represents the useful size of the NIR 
images. Table 2 also lists with which telescope the images 
were obtained, as well as the total on-source exposure time 
and the spatial resolution. The NIR images of all galaxies 
are shown in Fig. 1. 



4 BULGE/DISK DECOMPOSITION 

In recent years, performing bulge/disk (B/D) decomposi- 
tions in two dimensions, i.e., on the full image, instead of in 
one dimension, on the luminosity profiles, has become the 
accepted norm (e.g., Byun & Freeman 1995; de Jong 1996a; 
Mollenhoff & Heidt 2001; Simard et al. 2002; Tran et al. 



Galaxy 


Telescope 


Exp. time 


Seeing 


Size 






(min) 


(arcsec) 


(arcmin) 


NGC 0210 


WHT 


24 


1.04 


4.2 


NGC 0337A 


WHT 


28 


1.04 


4.2 


NGC 0488 


WHT 


45 


1.68 


4.9 


NGC 0628 


Bok 


N/A 


1.55 


11.0 


NGC 0864 


WHT 


24 


0.94 


4.2 


NGC 1042 


WHT 


45 


2.16 


4.8 


NGC 1068 


Bok 


N/A 


1.64 


9.6 


NGC 1073 


WHT 


24 


1.02 


4.2 


NGC 1169 


WHT 


24 


0.78 


4.2 


NGC 1179 


WHT 


28 


0.89 


4.2 


NGC 1300 


WHT 


32 


1.03 


5.2 


NGC 2775 


WHT 


12 


0.78 


4.2 


NGC 2805 


WHT 


48 


1.06 


4.2 


NGC 2985 


WHT 


20 


1.14 


4.2 


NGC 3184 


WHT 


24 


0.89 


5.1 


NGC 3227 


WHT 


12 


0.73 


4.7 


NGC 3344 


WHT 


24 


0.85 


4.9 


NGC 3351 


WHT 


12 


1.14 


7.4 


NGC 3368 


WHT 


8 


0.77 


5.0 


NGC 3486 


WHT 


24 


0.83 


4.9 


NGC 3631 


WHT 


38 


0.87 


4.2 


NGC 3726 


WHT 


12 


0.69 


4.2 


NGC 3810 


WHT 


12 


0.69 


4.2 


NGC 4030 


WHT 


30 


0.69 


4.2 


NGC 4051 


WHT 


12 


0.74 


4.8 


NGC 4123 


WHT 


28 


0.76 


4.2 


NGC 4145 


WHT 


36 


1.04 


4.2 


NGC 4151 


WHT 


12 


0.75 


4.8 


NGC 4242 


WHT 


24 


1.26 


4.2 


NGC 4254 


WHT 


12 


0.82 


4.2 


NGC 4303 


WHT 


19 


1.06 


4.8 


NGC 4314 


WHT 


35 


2.02 


4.2 


NGC 4321 


WHT 


12 


0.76 


5.0 


NGC 4395 


WHT 


20 


1.04 


4.8 


NGC 4450 


WHT 


12 


0.89 


4.7 


NGC 4487 


WHT 


13 


0.83 


4.2 


NGC 4535 


WHT 


12 


0.77 


5.0 


NGC 4548 


WHT 


12 


0.81 


4.8 


NGC 4579 


WHT 


12 


0.98 


4.8 


NGC 4618 


WHT 


28 


1.16 


4.2 


NGC 4689 


WHT 


12 


0.81 


4.2 


NGC 4725 


WHT 


24 


0.97 


6.0 


NGC 4736 


WHT 


19 


0.98 


5.2 


NGC 5247 


WHT 


12 


0.81 


4.8 


NGC 5248 


WHT 


12 


0.79 


4.8 


NGC 5334 


WHT 


16 


0.94 


4.2 


NGC 5371 


WHT 


12 


0.91 


4.2 


NGC 5457 


WHT 


6 


3.03 


4.2 


NGC 5474 


WHT 


16 


1.00 


4.2 


NGC 5850 


WHT 


12 


0.93 


4.2 


in Vjvy oyzi 


WTTT 

VV XI 1 


12 


1 02 


4 2 


NGC 5964 


WHT 


32 


0.82 


4.2 


NGC 6140 


Bok 


N/A 


1.74 


8.0 


NGC 6384 


WHT 


8 


0.70 


4.2 


NGC 6946 


Bok 


N/A 


2.31 


11.6 


NGC 7727 


WHT 


56 


0.82 


4.2 


NGC 7741 


WHT 


75 


0.90 


4.8 



Table 2. Properties of the final NIR images: telescope used to 
obtain the image (col. 2), total on-source exposure time in minutes 
(col. 3), seeing as measured from the final image (in arcsec; col. 4), 
and field of view on the final image (in arcmin, where images are 
assumed to be square; col. 5). 



© 2003 RAS, MNRAS 000, 1-9 



6 J. H. Knapen et al. 



Figure 1. K B images of all galaxies in the sample. North is up and East to the left in all images, and the scale of each individual image 
is indicated by the white bar in the bottom right corner, which represents a length of 1 arcmin. Indicative contours are shown at surface 
brightness levels of 17 and 15 mag arcsec -2 for all galaxies, except the following, for which contours are shown at 16 and 15 mag arcsec -2 : 
NGC 2775, NGC 3368, NGC 4030, NGC 4051, NGC 4254, NGC 4303, NGC 4535, NGC 4548, NGC 4736, NGC 5248, NGC 6384. 



2003). We have made 2D fits to our A'-band images using 
the GIM2D code (Simard et al. 2002). Unfortunately, due to 
limitations in the signal-to-noise ratio and image extent, and 
due to mosaicing imperfections, 2D fitting of our if-band 
images proved to be highly unstable. We therefore resorted 
to ID B/D fitting, pointing out that at least the resulting 
disk parameters (and most of the bulge parameters) will be 
equal within the errors of the fit in the ID and 2D cases (de 
Jong 1996a; Mac Arthur et al. 2003). In addition, for face- 
on galaxies ID and 2D fitting is nearly identical when no 
non-axisymmetric features (bars, arms) are fitted and error 
weighting is done consistently. 

We used the IRAF ellipse fitting tasks to fit radial pro- 
files to the surface brightness distribution in all images. We 
used sets of fixed position angle ellipticity for each galaxy, as 
taken from the RC3 and checked with our NIR and optical 
(Paper II) images. The resulting azimuthally averaged pro- 
files are shown in Fig. 2. The ID fitting technique we used 
is very similar to the one described in de Jong (1996a), with 
some modifications reflecting the investigations described by 
Mac Arthur et al. (2003). The Levenberg-Marquardt nonlin- 
ear least-squares fitting technique was used. For the disks we 
used exponential profiles, for the bulges Sersic profiles. Be- 
fore fitting, the model profiles were convolved with Gaussian 
seeing profiles with FWHMs equal to the FWHMs measured 
on foreground stars in the images. To reduce the effects of 
uncertainties in the seeing profiles, we discarded the inner 
few points, namely those with radii less than two tenths 
of the FWHM seeing value. Relative weights of the profile 
data points were set to reflect photon statistics, uncertain- 
ties due to sky measurement errors, seeing errors and errors 
due to unmodelled features. The latter correction is neces- 
sary because the light profiles contain features that are not 
modelled, like bars and spiral arms, which cause residuals 
in the fit (see e.g., de Jong 1996a, his fig. 3). Not taking the 
unmatched components into account will result in reduced 
X 2 values much larger than 1. We have used an uncertainty 
of 0.05 mag rms at each data point due to unmatched com- 
ponents in our weighting scheme. 

As shown by MacArthur et al. (2003), fitting Sersic pro- 
files with a free n value is very unstable and often results 
in non-physical results. We have used a similar grid search 
technique as them, fitting the profiles with fixed Sersic n 
values, with n ranging from 0.5 to 5 in steps of 0.1. Unlike 
MacArthur et al. (2003), we did not use a combination of 
an inner and a global \ 2 to select the best fit, but used only 
the x 2 °f the whole profile fit to select the best fit. 

All fits were repeated with the maximum sky error es- 
timate added and subtracted from the luminosity profile. 
This provides a reliable estimate of the uncertainties in the 
bulge and disk parameters, as it has shown before that sky 
errors dominate B/D decompositions (e.g., de Jong 1996a; 
MacArthur et al. 2003). The fits of many galaxies had large 
errors because of the limited field of view of the detector 
compared to the galaxy size and the limits of mosaicing in 
the near-IR. We therefore decided to exclude all disk param- 




20 40 60 

Radius (arcsec) 



Figure 2. Continued. 



eters if the uncertainty in their central surface brightness 
(Ho) or disk scale length (h) was larger than 30% or if the 
fitted scale length was larger than half the image size. Like- 
wise, bulge fits were rejected if the error in the effective ra- 
dius (-R c ,b) or in the surface brightness at this radius {fJ, e ,B) 
was larger than 30%, if 7? e ,B was smaller than the FWHM 
of the seeing, or if n < 0.8 (as is the case for NGC 3351 and 
NGC 7741 which have disturbed centres). This strict selec- 
tion leaves only 13 of the 57 galaxies with both reliable disk 
and bulge parameters. In Fig. 2 we show the fits for all those 
galaxies where either bulge or disk fit was reliable. 

Fitting results for those galaxies with at least a reliable 
bulge or a reliable disk fit are presented in Fig. 3 and Table 3. 
Fits are not plotted for those galaxies where both bulge and 
disk parameters could not reliably be determined. Bulge or 
disk parameters that were poorly determined are marked by 
a star in Fig. 2. 

In what follows, photometric parameters were corrected 
for Galactic foreground extinction following the precepts 
of Schlegel, Finkbeiner & Davis (1998). Distances to put 
measurements at an absolute scale were obtained from the 
Nearby Galaxies Catalog (Tully 1988). Disk surface bright- 
nesses were corrected to face-on values using the fully trans- 
parent correction of no = no — 2.51og(l — e), where the 
ellipticities e were taken from the RC3 but double-checked 
with the outer parts of the images. 

Many combinations between the various bulge and disk 
parameters can be investigated, but with so few reliable de- 
compositions we refer to previous work for most of them 
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Name 




,B 


Re,B 




n 


Mo 


h 




(mag arcscc ) 


(arcsec) 






(mag arcsec -2 ) 


(arcsec) 


NGC 0337A 


21. 


U °-0.08 


25.73_^9 





q+0.1 
■ y -o.i 


iq 70+0.30 
la. I »_q 32 


65.40+^-j 2 


NGC 0488 


16. 


oc+0.13 
JO _0.37 


— 2 .DO 


1 


o+O.l 
b -0.3 


16.64_°;« 


34.18±?;£ 


NGC 0628 


18 


94+0.03 
zzt _0.14 


22 44+0.80 
ZZ - 44 -3.20 


1 


c+0.1 

,O -0.1 


17.68lg;« 


82.42_5-86 


NGC 0864 


18. 


oe+0.34 
JD -1.09 


6.8it 6 : 3 7l 


3 


q+0.3 


17.47t°;°g 


31.86^j0 


NGC 1042 


18. 


64 +o.2i 




2 


1+0.2 
-0.4 


18.69l° ;? 




NGC 1169 


16 


42 +0.12 
^* z -0.13 


7 99+0.59 
1 - z -0.71 


1 


6+0.1 
• u -o.i 


16.81_oT 2 


33.49l 2 ;^ 


NGC 1179 


19. 


01 +0.07 
' ol -0.14 


5 54+0.34 
°- 04 -0.68 


1 


7+0.1 
' ' -0.1 


18.921°;°* 


31.57+ 2 ;^ 


NGC 1300 


16. 


c;c+0.03 
OD -0.11 


7 94+0.17 

1 - z -0.48 


1 


Q + 0.1 

,y -0.1 


18.051°;?° 


48.7ll 2 ; 3 7 « 


NGC 2775 


16. 


o 4 +0.26 
° -0.28 


9.27lJ-;_ 


2 


4+0.2 
-0.2 


15.95_° \l 


26.20_ 2 ' 17 


NGC 3227 


14. 


o 6 +0.ll 
°°-0.20 


2.86t°i 6 2 


1 


Q + 0.1 

• y -0.2 


16.27_°;^7 


25.6413-11 


NGC 3351 


15 


45+0.03 

^°-0.03 


9.46+^ 





O+0.1 

• s -o.i 


16 58+ 06 


45-44l 5 5 i° 


NGC 3631 


17. 


91+O.II 
-0.03 


10.30t°;°? 


1 


o+O.l 

• s -o.i 


18.65±g;J| 


64.581^1 


NGC 3810 


16. 


-,5+0.25 
lo -0.13 


3 40+O.59 


1 


5+0.2 


15.68tg;J 2 


15.96_J : ™ 


NGC 4051 


14 


otr + 0.18 
oo -0.19 


o 99+O.29 
J - z -0.34 


2 


rr+0.2 
,o -0.2 


17 24+0-1 5 

11 - z *-0.14 


48.2311^ 


NGC 4303 


15. 


n o+0.09 
uo -0.03 


4 46+0.19 
*- 4D -0.09 


1 


Q+0.1 

•°-o.i 


16 51+0- 06 


3q 54+3-72 


NGC 4535 


16 


1Q+0.21 
— 0.11 


3 go+0.42 


2 


1+0.2 

■ L -o.i 


17.57_°;°1 


63.45_H;£ 


NGC 4689 


19 


1Q +0.14 
-0.31 


14.79±_;_? 


2 


i+O.l 

-0.2 


17.78lg;S 


42.951*^ 


NGC 4725 


16. 


46+0.13 

^"-0.13 


i6.i6ij;l 


2 


k+0.1 

,o -0.1 


17 S0+°-l 4 


83.67_2i'.?2 


NGC 5247 


18 


OO+0.05 
88 -0.15 


no qo + 1-58 


2 


o+O.l 

•°-o.i 


18-22^30 


90 q7+l 6 - 96 


NGC 5371 


16. 


9n +003 
z -0.03 


6-o6±g : _; 


1 


O+O.l 

• 8 -o.i 


17 14+0.05 
l'- i4 -0.05 


51 44+ 6 - 52 
o± -^-10.73 


NGC 5457 


18. 


63+0.14 
-0.06 


21.82+ 2 -23 


2 


i+O.l 

■ L -o.i 


17 5n +0 06 

l'- OU -0.06 


128.07+i^ 


NGC 5474 


19. 


46+0.03 

w -0.03 


18.16±°,; 7 5 5 





Q + 0.1 

• y -o.i 


19 7 i+0.l3 
ly - (1 -0.13 


8i-29l 2 30 :S 


NGC 5850 


17. 


cn+0.14 
Ou -0.25 




2 


q+0.1 

,y -0.2 


18.25±°;^ 


39.1lll°;_| 


NGC 5921 


15 


54 +0.12 
J -0.22 


3 5O+0-27 
J - oy -0.50 


1 


g+o.i 


1 7 O7+0.23 
1 1 ..31 _ 26 


30.43__23.89 


NGC 6140 


19. 


7O+0.06 
' -0.08 


, _i_4 fifi 

99 74i- 4 -oo 
' -6.34 


1 


i +0.1 

■ 1 -o.i 


. . in Ad 

1 q 03+ 


43.45+17°^ 


NGC 6384 


16 


01 +0.27 
Si -0.15 


q qq+l- 83 


1 


q+0.2 

• y -o.i 


16 96+ - 24 

10.30_ 20 


41 09+ 9 ' 27 


NGC 6946 


16 


qq+0.22 
uy -0.04 


12.24^ ;« 


2 


9+0.2 
-0.1 


16.95+^ 


127.18l 21 2 5° 


NGC 7727 


16. 


,o+0.14 
18 -0.24 




2 


4+0-1 

-0.2 


17.28t°;» 


29.581^2 


NGC 7741 


18 


oo+0.04 
88 -0.04 


8-43±8:ii 





g+o.i 


18.38_g;S? 


38.76+Hl 



Table 3. Fitted bulge and disk parameters, not corrected for extinction or inclination 



(e.g., Moriondo ct al. 1998; Mollenhoff & Heidt 2001; Gra- 
ham 2001). We note, however, that for spiral galaxies cor- 
relations between the disk parameters are severely affected 
by selection effects, most notably in the case of correlations 
involving ^io and h. Intrinsically large, high surface bright- 
ness galaxies are simply more likely to enter a magnitude 
or diameter limited sample than small, low surface bright- 
ness galaxies (Disney & Phillipps 1983; Allen & Shu 1979). 
These selection effects have to be taken into account to cal- 
culate proper distributions (de Jong & Lacey 2000; Blanton 
et al. 2001; Cross & Driver 2002). Furthermore, if bulge pa- 
rameters are tightly correlated with disk parameters (like 
the h/ 7? e ,B ratio discussed below), then also the distribu- 
tion of bulge parameters will be biased. 

In Fig. 3 we show the reliable bulge and disk param- 
eters as function of morphological T-type (from the RC3) 
and of best fitting n. As noted before (e.g., de Jong 1996b; 
Mollenhoff & Heidt 2001), the change from T-type 2 to 8 
is strongly correlated with the change in bulge luminosity. 
The decrease in luminosity with increasing T-type is mainly 



the result of bulge surface brightness changes helped a little 
bit by bulge shape changes, but is not due to bulge scale 
size changes. 

Figure 3 also shows that the disk luminosities peak at 
about T-type 3-4, decreasing on average for earlier and 
later-types, driven mainly by surface brightness but partly 
by scale size, as observed before by de Jong (1996b) and 
Mollenhoff & Heidt (2001). The relative bulge to disk scale 
size ratio does not change much in this morphological 
type range as often observed before (e.g., de Jong 1996b; 
Courteau et al. 1996; Moriondo et al. 1998; Graham 2001; 
Mac Arthur et al. 2003). This implies that the change ob- 
served in B/T ratio as a function of morphological type, 
a change which is in itself not surprising as the B/T ra- 
tio is one of the criteria for defining morphological type, 
is mainly due to differences in the relative effective surface 
brightnesses of the bulge and the disk. 

When now turning to the bulge and disk parameter cor- 
relations with n, we confirm previously known trends be- 
tween n and bulge luminosity and surface brightness (e.g., 
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Parameters 


r s 


Significance 


Type - B/T 


-0.75 


0.047 


Type - ^t c>B 


0.46 


0.169 


Type - M B 


0.81 


0.002 


n - M B 


-0.54 


0.083 


Type - n 


-0.53 


0.090 



Table 4. For a number of selected correlations, as discussed in 
the text (column 1), we list the Spearman rank-order correlation 
coefficient r s (column 2) and significance (column 3), i.e., the 
probability that this relation occurs at random. 

Andredakis et al 1995; Mollenhoff & Heidt 2001; Trujillo et 
al. 2002). These authors also report other correlations be- 
tween n and, for instance, R c ,b and B/T, but we have too 
few data points to confirm such trends. Finally, there is a 
clear trend between morphological type and n, observed by 
most other authors, though not as strongly by Mac Arthur 
et al. (2003). We present correlation coefficients and signifi- 
cance indicators of some of the trends mentioned in Table 4. 
Almost all trends and parameter ranges are consistent with 
continuations of trends seen in a large sample of early-type 
galaxies studied by R. S. de Jong et al. (in preparation). 

5 CONCLUDING REMARKS 

In this paper, we present the first set of results of our study of 
the spiral arm and disk properties of a sample of 57 nearby, 
relatively face-on, spiral galaxies. We describe the sample se- 
lection and summarise some of the properties of the sample, 
and present a full set of NIR Jf s -band images of the sample 
galaxies. Those images have been obtained using the 4.2 m 
WHT and 2.3 m Bok telescopes, and with a minimum im- 
age size of 4.2 arcmin square cover most of the disk of all 
galaxies. We describe the data collection and reduction pro- 
cedures, and show all resulting images. 

We derive radial profiles to the light distribution of all 
sample galaxies, and perform ID bulge/disk decomposition 
by fitting exponential disk and Sersic bulge profiles. The 
results of this exercise are shown, but due to the limited 
signal-to-noise ratio in our images we can fit both bulge and 
disk parameters in a reliable way for only 13 of our 57 galax- 
ies, although either reliable bulge or disk parameters can be 
determined for more galaxies. Using the reliable bulge and 
disk parameters thus derived, we explore several correlations 
between those and galactic parameters such as morpholog- 
ical type. We confirm correlations between bulge and disk 
luminosities and morphological type, and between the Ser- 
sic parameter n and bulge luminosity and surface brightness, 
and morphological type. 

In further papers in this series, we will publish the 
accompanying optical broad (B, I) and narrow-band (Ha) 
imaging, as well as analyses of disk scale lengths, and of 
properties of morphological entities such as bars and rings 
within the galaxies. 
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